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The crystal structure and local spontaneous polarization of (BaTi03) m / (SrTiOs)™ superlattices is 
calculated using a first-principles density functional theory method. The in-plane lattice constant is 
1% larger than the SrTi03 substrate to imitate the relaxed superlattice structure and the symmetry 
is lowered to monoclinic space group Cm which allows polarization to develop along the [110] and 
[001] directions. The polarization component in the [110] direction is found to develop only in the 
SrTi03 layers and falls to zero in the BaTiC>3 layers, whereas the polarization in the [001] direction 
is approximately uniform throughout the superlattice. These findings are consistent with recent 
experimental data and first-principles results for epitaxially strained BT and ST. 

PACS numbers: 68.03.Hj, 68.65.Cd, 77.80.-c, 81.05.Zx 



I. INTRODUCTION 

Perovskite materials, such as BaTiC>3 (BT) and 
SrTiC>3 (ST), arc of considerable interest due to their 
polarization-related properties. Short-period superlat- 
tices of these materials have switchablc polarization, fer- 
roelectric T c and dielectric and piezoelectric coefficients 
distinct from those of their individual constituents, and 
the freedom to choose from a variety of constituent ma- 
terials, layer thicknesses and sequences can be exploited 
to optimize selected desirable properties 0, 0, IE El IE IE 
ill. 

The (BaTi0 3 )m-(SrTi0 3 )„ (BT/ST m/n) superlattice 
has been extensively investigated by experimentalists pri- 
marily because of the possible enhancement of the switch- 
able polarization arising from strain effects. Due to the 
lattice mismatch with ST, there is a compressive in-plane 
strain on BT that increases the tetragonal strain and, 
through polarization-strain coupling, the [001] polariza- 
tion. Several groups Q, |E El studied the structure 
and dielectric properties of symmetric and asymmetric 
BT/ST superlattices with periods ranging from 4 to 500 
unit cells. Measurements of these superlattice structures 
were consistent with tetragonal symmetry. Depending 
on the BT/ST ratio, the polarization can be enhanced 
over that of unstrained BT and the Curie temperature 
elevated. 

Neaton and Rabe Q performed first-principles calcu- 
lations to study the structure and polarization of BT/ST 
superlattices lattice-matched to an ST substrate. With 
this epitaxial constraint the systems studied (m/n su- 
perlattices with m + n = 5) were found to be stable with 
PAmm symmetry, with polarization along [001]. The to- 
tal polarization was decomposed into contributions from 
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individual Ti-centred layers. It was found that the un- 
strained ST layers were polarized along [001] and the 
strain-enhanced polarization of BT reduced such that the 
polarization was nearly constant in the [001] direction. 
This was attributed to the minimisation of electrostatic 
energy costs associated with the buildup of polarization 
charge at the interfaces, suggesting a simple model that 
provided a good fit to the first-principles results. 

Recently Jiang et al. || and Rfos et al. ^l| investi- 
gated the structure of BT/ST superlattices with periods 
ranging from 4 to 100. Through a combination of second 
harmonic generation and x-ray diffraction measurements, 
they found that below a critical periodicity of 30/30, the 
system has symmetry lower than tetragonal, consistent 
with a nonzero polarization in the ST layers along the 
[110] direction. The in-plane lattice constants are «0.8% 
larger than that of bulk ST which means that the films 
are at least partially relaxed. 

In this paper we describe first-principles calculations 
that explore the idea that it is precisely this increase in 
in-plane lattice constant that is responsible for the ob- 
served change in structure. For several BT/ST superlat- 
tices with period up to 10, we expand the in-plane lattice 
constant to 1.01 times the theoretical lattice constant of 
bulk ST and relax the tetragonal symmetry constraint 
to allow polarization to develop along the [110] direction 
as well as along [001]. We show that the ST layers ex- 
hibit both the [001] polarization which, as in the coherent 
tetragonal superlattices, minimises the electrostatic en- 
ergy due to charge accumulation at the BT-ST interface 
layers, as well as a nonzero component of the polarization 
along [110], lowering the symmetry as observed. The de- 
tails of the methods used are given in Section [H] and in 
Section ITTT1 we present our results and discussion. 



II. METHOD 

First-principles calculations were performed with the 
Vienna ab initio Simulations Package (VASP) ^2] which 
implements density functional theory within the local 
density approximation. Projector augmented wave po- 
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tentials were used ^301 with the semicore Ba (5s, 5p), 
Sr (4s,4p), Ti (3s, 3p) treated as valence states. A plane 
wave energy cutoff of 600 eV and a 6 x 6 x 2 Monkhorst- 
Pack fc-point sampling of the Brillouin zone were used. 

The super lattices have structure m/n where m and n 
represent the thickness in unit cells of BT and ST respec- 
tively. The structures considered in this paper are 3/3, 
4/4. 4/6, 5/5 and 6/4, with unit cells 1 x 1 x (n + rri). 
The structures are optimized in the monoclinic space 
group Cm, which allows nonzero polarization of the form 
P\(x + y) + P3S. Because the deviation from the in- 
plane square lattice observed in the x-ray experiment is 
very small[Tl|. we constrain the angles between the lat- 
tice vectors to 90 degrees. The c lattice parameter and 
the ions are allowed to fully relax within the symmetry 
of the Cm space group. The in-plane lattice parame- 
ter is fixed to a — I.OIost = 3.902 A, with the calcu- 
lated values of the lattice constants for the constituent 
materials being agT = 3.863 A for bulk cubic ST and 
a B T = 3.943 A and c/a = 1.012 for bulk tetragonal BT. 
This results in an epitaxial strain of -1.04% on the tetrag- 
onal BT layer. The relaxation was considered to be com- 
plete when the Hellman-Feynman forces on the ions were 
less than 5 meVA -1 . 

As discussed in the Introduction, this value of in-plane 
lattice constant was chosen to mimic the experimental re- 
laxation of the superlattice. This relaxation implies that 
the BT/ST superlattices, although epitaxial, are not co- 
herent with the substrate and there must be dislocations 
present. It is assumed that these dislocations are confined 
to the region near the substrate interface and therefore 
we are modelling the relaxed, defect-free superlattice far 
from the substrate interface. 

To investigate the spatial variation of the polarization 
and strain within the superlattice, we compute the lo- 
cal polarization in each 5-atom cell of the superlattices, 
defined as 

-Plocal = ^{Z^Az Ti + Z^Azoi + Zq 2 Az 2 

+ ±(Z* 03 _Az 03 - + Z* 03+ Az 03+ )} 

where Z* is the dynamical charge of atom i [2l| and Az; is 
the displacement in the z direction of atom i with respect 
to the centrosymmetric cell defined by the Ba/Sr atoms. 
03+ and 03_ represent the 03 atoms at the top and 
bottom of the 5-atom cell, respectively. The local strain 
c/a is obtained by using the A atoms to define the corners 
of the local unit cell, i.e. c/a = [z(A+) — z(A_)] /a. 

ILL RESULTS AND DISCUSSION 

The local c/a ratio for each layer is given in Table [I] 
The BT layers, being under compressive in-plane strain, 
have c/a larger than the bulk tetragonal value, whereas 
the ST layers are under tensile in-plane strain and have 
c/a < 1. The two interface layers are inequivalcnt, with 



TABLE I: Variation of c/a in each superlattice. The average 
c/a ratios for the BT and ST layers are shown at the bottom 
of the table. 
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the polarisation being directed towards the Ba in one 
and towards the Sr in the other. Both interfaces have 
c/a values intermediate between the average BT and ST 
values. 

Figures ^and |21 show the variation of Pi OC ai in each su- 
perlattice. The [001] component remains approximately 
uniform throughout the superlattice, dipping slightly in 
the ST layers. As the relative thickness of the BT layer 
increases, the magnitude of the [001] polarization in- 
creases. In the ST layers, there is a nonzero polariza- 
tion along [110], lowering the symmetry from PAmm. It 
grows smoothly from zero at the interface over about two 
lattice constants, reaching its maximum value in the cen- 
ter of the layer. As the relative thickness of the ST layer 
increases, the magnitude of the polarization along [110] 
increases. The [001] and [110] components are not in- 
dependent. This can be seen by comparing the 4/4 and 
6/4 superlattices, in which the ST layer has the same 
thickness. Relative to the 4/4 superlattice, in the 6/4 
superlattice the [001] component is larger and the [110] 
component is smaller. An analogous behavior is seen in 
the 4/4 and 4/6 superlattices. In the 4/6 superlattice, 
the [001] component is smaller and the [110] component 
is larger. In all superlattices studied except the shortest 
period (3/3), the magnitude of the total polarization in 
the center of the ST layer is quite similar. 

The local variation of polarization and strain can be 
understood by applying the model discussed in Ref. to 
the relaxed superlattice. The first step is to consider the 
effect of epitaxial strain on the polarization and c/a. In 
the relaxed BT/ST superlattice the BT layers are sub- 
ject to compressive in-plane strain whereas the ST lay- 
ers are subject to tensile in-plane strain. Dieguez et al. 
|15| studied the epitaxial strain phase diagram of bulk 
BT in zero macroscopic field using a first-principles ap- 
proach. They found that for compressive strains greater 
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FIG. 1: (Color online) Magnitude and direction of Pi oca i in 
the various superlattices. 
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FIG. 2: (Color online) [001] and [110] components and mag 

nitude Of Plocal • 



than 0.64% BT is in the c phase (tetragonal structure 
P4mm), in which the polarization points along the [001] 
direction. By extrapolating their structural parameters 
to the relaxed lattice constant a=3.902 A we calculate the 
local polarization and local c/a of BT to be 28.4^Ccm -2 
and 1.036, respectively. A similar study of the effect 
of strain on ST was carried out by Antons et al. |l6j |. 
They found that for tensile strains greater than 0.54% 
ST has orthorhombic structure Amm2 with polarization 
in the [110] direction. The predicted polarization and 



c/a of bulk ST with an in-plane tensile strain of 1% are 
« 20/xCcm~ 2 and 0.985, respectively [Tt| . 

The next step is to take into account the electrostatic 
energy of combining these two strained layers into a su- 
perlattice. In particular, a discontinuity in the [001] com- 
ponent of the polarization is energetically costly. The en- 
ergy is reduced by concomitant polarization of the highly 
polarizable ST layer and reduction in the polarization of 
the BT layer to achieve continuity. The final value of the 
polarization is strongly dependent on the relative thick- 
nesses of the two layers. The thicker the BT layer rela- 
tive to the ST layer, the higher the polarization. These 
trends are evident in the computed polarizations shown 
in Figure |5J 

Through polarization-strain coupling, these changes in 
the polarization lead to changes in the c/a ratio. This 
can at least roughly account for the differences of the 
local BT and ST c/a from the strained bulk values, and 
for the trends between c/a and polarization in Table 

In the ST layer, increasing polarization along [001] re- 
duces the component along [110], minimizing changes in 
the total magnitude. This coup ling has previously been 
noted in a first-principles study |l6l ] of phonon modes in 
strained ST, where in the Amm2 phase, nonzero polar- 
ization along [110] is associated with a hardening of the 
ferroelectric mode along [001]. This observation allows 
us to account for the observed variation in the maximum 
[110] component in Figure |2 

In comparing these results with experiment, it is im- 
portant to select only systems with in-plane relaxation, 
such as the superlattices in Ref. 0, 0. The shortest pe- 
riod system studied there is the 10/10 superlattice. We 
can extrapolate the results from the three symmetrical 
superlattices, namely the 3/3, 4/4 and 5/5. The c/a val- 
ues for each superlattice are remarkably similar so we 
can expect the 10/10 superlattice to have similar c/a 
values, as indeed it does. We can estimate the period 
of the 10/10 superlattice by scaling the periods of the 
3/3, 4/4 and 5/5 superlattices. This gives « 78.5 A 
which compares well with measured values of 78.1 A [ll| 
and 79 A j6j. Consistent with the simple electrostatic 
model, the [001] component increases only very slightly 
with increasing period, and saturates to a estimated value 
of w 22 — 23 /iCcm 2 , about 90% percent of the com- 
puted free bulk polarization of BT The [110] com- 
ponent at the centre of the ST layers increases from less 
than 10 /iCcm" 2 for the 3/3 to « 13^Ccm" 2 for 4/4 to 
f=a 15^tCcm~ 2 for the 5/5. This should be close to the sat- 
uration value; smaller than the strained bulk ST polar- 
ization of w 20/iCcm~ 2 , as expected, due to the coupling 
to the nonzero [001] component. 

To extrapolate to still longer periods, it is necessary to 
recognize that the behavior will change once the thickness 
of the layers exceeds the polarization coherence length 
and/or the elastic critical thickness of the individual lay- 
ers. The uniformity of the [001] polarization is expected 
only for short period superlattices, as for thicker layers 
the energy cost of nonzero div P competes with the en- 
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ergy cost for polarizing the paraelectric ST. As a result, 
the [001] component decays into the interior of the ST 
layers and for sufficiently thick layers will become zero. 
This behavior is demonstrated in KNbOs/KTaO^ super- 
lattices shell model studies by Sepliarsky et al. El ■ In 
addition, as the individual layers exceed their elastic crit- 
ical thickness, they would be expected to relax to their 
separate in-plane lattice constants through the formation 
of misfit dislocations. The precise value depends on the 
method of synthesis; it has been shown that coherent 
films can be grown well over the critical thickness under 
appropriate conditions 0, 0] . 

In summary it seems that the overall strain state is 
crucial to the properties of the superlattice. If the su- 
perlattice relaxes, as in the work of Rios et al. the ferro- 



electric properties can be dramatically different than in 
a superlattice where the strain is preserved. For BT/ST, 
a simple model of electrostatics and bulk strain effects 
works very well, with the very gentle interface between 
BT and ST being relatively unimportant. 



Acknowledgments 

We would like to thank Refik Kortan, Susana Rios, 
Andreas Rudiger and Jim Scott for valuable discussions. 
This work was supported by DOE Grant DE-FG02- 
01ER45937 and ONR Grant N00014-00-1-0261. 



[1] H. Tabata, H. Tanaka, and T. Kawai, Appl. Phys. Lett. 

65, 1970 (1994). 
[2] M. Sepliarsky et al, Phys. Rev. B 64, 060101 (2001). 
[3] M. Sepliarsky et al, Journal of Applied Physics 90, 4509 

(2001) . 

[4] H. M. Christen, E. D. Specht, S. S. Silliman, and K. 

Harshavardhan, Phys. Rev. B 68, 020101 (R) (2003). 
[5] A. Jiang, J. Scott, H. Lu, and Z. Chen, J. Appl. Phys. 

93, 1180 (2003). 
[6] T. Shimuta et al, Journal of Applied Physics 91, 2290 

(2002) . 

[7] J. Kim et al., Integrated Ferroelectrics 47, 235 (2002). 
[8] J. Neaton and K. Rabe, Appl. Phys. Lett. 82, 1586 

(2003) . 

[9] H. Tabata and T. Kawai, Appl. Phys. Lett. 70, 321 
(1997). 

[10] H. Tabata, H. Tanaka, T. Kawai, and M. Okuyama, Jpn. 

J. Appl. Phys. I 34, 544 (1995). 
[11] S. Ri'os et al, J. Phys.: Condens. Matter 15, L305 (2003), 

bST. 

[12] G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 
(1996). 



[13] P. Blochl, Phys. Rev. B 50, 17953 (1994). 

[14] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

[15] O. Dieguez et al., Phys. Rev. B 69, 212101 (2004). 

[16] A. Antons, J. Neaton, K. M. Rabe, and D. Vanderbilt, 
http: / /arxiv.org/abs /cond-mat /0407077 (2004) . 

[17] J. Haeni et al., Nature 430, 758 (2004). 

[18] K. Choi et al., Science 306, 1005 (2004). 

[19] P. Ghosez, X. Gonze, and J. -P. Michenaud, Ferroelectrics 
206-207, 205 (1998). 

[20] W. Zhong, R. King-Smith, and D. Vanderbilt, Phys. Rev. 
Lett. 72, 3618 (1994). 

[21] The dynamical charges are approximated by those of bulk 
cubic BT in the BT layers, bulk cubic ST in the ST 
layers and average values in the two Ti02 interface layers. 
The values of the dynamical charges for cubic ST and 
tetragonal BT were taken from Refs. [Tills The variation 
of -Piocai depends on the choice of Z* . Even though the 
effect is small it should be noted that Pi OC ai is only a 
guide to the variation of the polarization and should not 
be treated as an exact quantity. 



